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(57) Abstract: Microwave imaging via space-time beamforming is carried out by transmitting microwave signals from multiple 
antenna locations into an individual to be examined and receiving the backscatteied microwave signals at multiple antenna locations 
to provide received signals from the antennas. The received signals are processed in a computer to remove the skin interface reflection 
component of the signal at each antenna to provide corrected signal data. The corrected signal data is provided to a beamformer 
process that time shifts the received signals to align the returns from a scatterer at a candidate location, and then passes the time 
aligned signals through a bank of filters^ the outputs of which aie summed, time-gated and the power therein calculated to produce 
the beamformer output signal at a candidate location. The beamformer is then scanned to a plurality of different locations in the 
individual by changing the lime shifts, filter weights and time-gating of the beamformer process. The output power may be displs^red 
as a function of scan location, with regions of lai^e output power corresponding to significant microwave scatterers such as malignant 
lesions. 
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SPACE-TIME MICROWAVE IMAGING FOR CANCER 

DETECTION 

FIELD OF THE INVENTION 

[0001] The present invention pertains generally to the field of medioal 
imaging and particularly to microwave imaging of tissue for the detectiori 
and location of tumors. 

BACKGROUND OF THE INVENTION 

[0002] Various imaging techniques have been employed for detecting 
and locating cancerous tumors in body tissue. X-ray and ultrasound 
imaging techniques are commonly utilized in screening for breast cancer. 
X-ray mammography is the most effective current method for detecting 
early stage breast cancer. However, X-ray mammography suffers from 
relatively high false positive and false negative rates, requires painful 
breast compression, and exposes the patient to iow levels of ionizing 
radiation. 

[0003] Microwave based imaging methods have been proposed for use 
in imaging of breast tissue and other body tissues as an alternative to 
current ultrasound and X-ray imaging techniques. Microwave imaging 
does not require breast compression, does not expose the patient to 
ionizing radiation, and can be applied at low power levels. Microwave- 
based imaging exploits the large contrast in dielectric properties between 
normal and malignant tissue. With microwave tomography, the dielectric 
properties profile of an object being imaged is recovered from 

■ 

measurement of the transmission of microwave energy through the 
object. This approach requires the solution of an ill-conditioned nonlinear 
inverse-scattering problem which requires elaborate image reconstruction 
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algorithms. An alternative microwave imaging approach is based on 
backscatter methods that use the measured reflected signal to Infer the 
locations of significant sources of scattering in the object being imaged, 
and are simpler to implement and more robust. Backscatter methods 
require the focusing of the received signal in both space and time to 
discriminate against clutter and to obtain acceptable resolution. This may 
be accomplished with an antenna array and ultra-wideband microwave 
probe signals. For a discussion of this approach, see, S.C. Hagness, et 
aL, "Two-Dimensional FDTD Analysis of a Pulsed Microwave Confocal 
System for Breast Cancer Detection: Fixed Focus and Antenna-Array 
Sensors," IEEE Trans. Biomed. Eng., Vol. 45, Dec, 1998, pp. 1470- 
1479; S-C. Hagness, et aL, "Three-Dimensional FDTD Analysis of a 
Pulsed Microwave Confocal System for Breast Cancer Detection: Design 
of an Antenna-Array Element," IEEE Trans. Antennas and Propagation, 
Vol. 47, May, 1999, pp. 783-791; S.C. Hagness, et aL, "Dielectric 
Characterization of Human Breast Tissue and Breast Cancer Detection 
Algorithms for Confocal Microwave Imaging," Proc. of the 2"** World 
Congress on Microwave and Radio Frequency Processing, Orlando, 
Florida, April, 2000; and X. Li, et aL, "A Confocal Microwave Imaging 
Algorithm for Breast Cancer Detection," IEEE Microwave and Wireless 
Components Letters, Vol. 11, No. 3, March, 2001, pp. 130-132. 

SUWIWIARY OF THE INVENTION 

[0004] Space^ime microwave imaging in accordance with the invention 
overcomes many of the limitations of conventional breast cancer 
screening modalities. Of particular significance, microwave imaging via 
space-time (MIST) beamforming has the potential to reduce the number of 
false negatives currently associated with conventional X-ray 
mammography. MIST exploits the significant dielectric-properties 
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contrast between malignant and normal breast tissue at microwave 
frequencies by taking advantage of the exceptionally strong biophysical 
contrast mechanisms of clinical interest, such as water content, 
vascularization/angiogenesis, blood flow rate, and temperature. MIST 
techniques in accordance with the invention have the potential for 
sensitivity and resolution sufficient to allow reliable detection of extremely 
srnall {millimeter size) malignant tumors even in radiographically dense 
breast tissue or in the upper outer breast quadrant near the chest wall- In 
addition, MIST utilizes non-ionizing microwave radiation, is noninvasive, 
does not require the injection of contrast agents, and avoids the need for 
breast compression. The present invention also has the potential to 
reduce the number of false positives associated with conventional X-ray 
mammography and thereby reduce the number of unnecessary biopsies. 
Because low-power microwave exposure is harmless, MIST exams may 
be done more frequently than X-ray mammography, and monitoring and 
comparison of breast scans from one exam to the next can be used to 
identify changes in lesions due to vascularization and the growth of 
cancerous tissue. Further, discrimination between malignant and benign 
tumors may also be possible based on spectral and polarization 
characteristics of benign and malignant tumors, MIST may be 
implemented utilizing relatively low-cost hardware, allowing reduced cost, 
screening procedures and allowing routine screening to be made more 
widely available to medically under-served populations in both developed 
and underdeveloped countries. Further, the safety of MIST imaging 
techniques, the comfort of the procedure (no breast compression 
required), the ease of use, and the low cost of the scanning procedure 
should help to improve acceptance by the public of regular (e.g., annual) 
screenings. The present invention further provides significantly improved 
imaging over prior microwave imaging approaches using a single antenna 
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with no spatial focusing or using an array of antennas with simple spatial 
focusing via time shift and summing of backscattered waveforms. 

tO0O5] The space-time microwave imaging system of the invention 
preferably carries out processing to remove the effect of artifacts such as 
the skin response of the antenna response from the data before beam 
former processing. The artifact removal process estimates the artifact 
component of the signal at each antenna in an array of antennas as a 
filtered combination of the signals at all other antennas. The filter 
weights are chosen to minimize the residual signal over that portion of the 
received data dominated by the artifact* 

[0006] The beam forming processing preferably first time shifts the 
received signals from the antennas after artifact removal to approximately 
align the returns from a scatterer at a candidate location. The time- 
aligned signals are passed through a bank of filters {e.g., finite-impulse 
response filters for time domain processing), one for each antenna 
channel, with the outputs of the filters then summed and time gated and 
the power calculated to produce the beamf ormer output signal at a 
candidate location. The filters may be designed using a least squares 
technique to present maximum gain (e.g., unit gain) to scattered signals 
originating from the candidate location. The beamf ormer output power 
thus represents an estimate of the energy scattered by that location. The 
beamformer process is then scanned to a plurality of different locations in 
the individual (e.g., in the breast) by changing the time shifts, filter 
weights, and time gating in the beamformer process. The output power 
may then be displayed as a function of scan location, with regions of 
large output power corresponding to significant microwave scatterers 
(e.g., malignant lesions). Processing may further be carried out to 
account for the effect of frequency dependent scattering. Scattering is 
frequency dependent due to dispersive dielectric properties and the 
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presence of multiple scattering surfaces. Errors due to frequency 
dependent scattering may be compensated by processing the beamformer 
output signal prior to time gating using parametric signal processing 
models for frequency dependent scatting effects. 

[0007] The present invention may also be utilized to carry out 
hyperthermia treatment of a detected lesion by applying signals to the 
antennas in the array with appropriate weights focus the microwave 
radiation from the antennas onto the detected position of the lesion. 

[0008] A space-time microwave imaging system that carries out 
imaging in accordance with the invention includes an array of antennas 
for radiating and receiving microwaves, a microwave source connected to 
the array of antennas to provide microwave signals such as pulse signals 
of a selected width and repetition rate to the antennas, and a receiver 
connected to the antennas to detect the microwave signals received by 
the antennas and provide signal data corresponding thereto. The system 
of the invention may also utilize a microwave source which provides 
discrete frequency signals that can be combined to provide the effect of a 
broadband pulse source. A computer is connected to receive the signal 
data and to carry out beamformer processing. The computer is also 
preferably programmed to estimate an artifact reflection component of a 
signal at each antenna as a filtered combination of the signals at a|l other 
antennas and to subtract the estimated artifact reflection component from 
the signal data to provide corrected signal data. The weights of the filters 
are chosen to miirimize a residual signal over that portion of the received 
data dominated by the reflection. The computer is programmed to 
process the corrected signal data in a beamformer process to time shift 
the corrected signal data to approximately align the returns from a 
scatterer at a candidate location. The time aligned signals are passed 
through a bank of filters with a filter for each antenna, the outputs of the 
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filters are summed to form a summed signal, and the power in the 
summed signal is calculated to produce a beamformer output signal. The 
beamformer filters are designed to present maximum gain to scattered 
signals originating from the candidate location. The beamformer process 
is scanned to a plurality of different candidate locations in the object to be 
imaged, such as a portion of an individual, by changing the time shifts 
and filter weights to generate multi-dimensional output data. An output 
device such as a cathode ray tube, LCD screen, etc. may be connected to 
the computer to display the multi-dimensional output power as a function 
of scanned locations^ providing an image on which cancerous lesions may 
be distinguished from surrounding tissue. The computer may be further 
programmed to time gate the summed signal to form a time-gated 
summed signal and to calculate the power in the time-gated summed 
signal. The beamformer process can be scanned in the object by 
changing the time shifts, filter weights and time gates to generate the 
multi-dimensional output power data. The computer may be further 
programmed to process the beamformer output signal from the filters 
prior to time gating using a parametric signal processing model to 
compensate frequency dependent scattering effects. The computer may 
be further programmed in the beamformer process to apply a selected 
window to the time aligned signals before passing the time aligned signals 
through a bank of filters, and to apply a selected window to the summed 
signal before the power in the summed signal is calculated, to reduce the 
effects of clutter in the signal. The beamformer filters are preferably FIR 
filters designed to satisfy a penalized least squares condition to present 
unit gain to scattered signals originating from a candidate location. To 
increase the signal to clutter ratio, the system may include signal ^ 
processing circuitry that receives the pulses from the microwave source 
and passes the pulses through a delay and a filter for each antenna before 
providing the delayed and filtered puises to the antennas/ The .delays and 
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filters for each antenna are selected to focus the radiated microwave 
energy from the array of antennas at a selected candidate location in the 
object. At sufficiently high power levels, such focussing may also be 
utilized to provide hyperthermia treatment at a location which has been 
previously identified as a tumor. 

[0009] Further objects, features and advantages of the invention will be 
apparent from the following detailed description when taken in 
conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0010] In the drawings: 

[0011] Fig- 1 is a block diagram of a space-time microwave imaging 
system in accordance with the invention for transmitting and receiving 
using the same antenna. 

[0012] Fig, 2 is a block diagram of a further embodiment of a space- 
time microwave imaging system in accordance with the invention 
providing simultaneous transmission and reception with all antennas. 

[0013] Fig. 3 is an illustrative view of an antenna array and its 
utilization in the microwave imaging system of the invention. 

[0014] Fig. 4 are spatial waveforms illustrating pulses reflected from 
the skin. 

[0015] Fig. 5 is an illustrative temporal waveform at an antenna 
showing reflected pulses from the skin interfaces. 

[0016] Fig, 6 is a block diagram illustrating the process of artifact 
removal from a backscattered signal at the first antenna (antenna 1). 
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[0017] Fig, 7 is a block diagram of a space-time beamformer for a 
location utilizing time-domain processing. 

[0018] Fig. 8 are time waveforms showing skin artifact removal, with 
the solid curve showing the original waveforms and the dashed curves 
corresponding to the waveforms after application of the skin artifact 
removal algorithm. 

[0019] Fig. 9 is a block diagram illustrating a preferred process for 
accommodating frequency-dependent scattering. 

[0020] Fig. 10 is a block diagram of another space-time beamformer for 
a location r^ utilizing frequency domain processing. 

[0021] Fig. 1 1 is a block diagram as shown in Fig. 6 with an addition of 
a component to reduce distortions from the skin response removal 
process. 

[0022] Fig. 12 is a block diagram for transmission of microwave pulses 
from an array of antennas to focus the microwaves at a chosen location, 
and which may be utilized In hyperthermia treatment of tumors. 

. DETAILED DESCRIPTION OF THE INVENTION 

[0023] In one embodiment for carrying out space-time microwave 
imaging (MIST) in accordance with the invention, each antenna in an 
array of antennas sequentially transmits a low-power ultra-short 
microwave pulse into an object to be imaged, such as the breast, and 
collects the backscatter signal. The relative arrival times and amplitudes 
of backscattered signals received by the antennas across the antenna 
array provide information that can be used to detect the presence and 
determine the location of malignant lesions. Breast carcinomas act as 
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significant microwave scatterers because of the large dieiectric-^properties 
contrast with the surrounding tissue. The problem of detecting and 
localizing scattering objects using pulsed signals and antenna arrays is 
similar to that encountered in radar systerns^ such as those used for air 
traffic control, military surveillance, and land-mine detection. 

[0024] Data in published literature and from our measurements on 
freshly excised breast biopsy tissue suggest that the malignant-to-normal 
breast tissue contrast in dielectric constant, s^, and conductivity, a, is 

between 2:1 and 10:1, depending on the density of the normal tissue. 
The higher dielectric properties of malignant breast tissue arise, in part, 
from increased protein hydration and a breakdown of cell membranes due 
to necrosis. The contrast ratio does not vary significantly with tumor age, 
which suggests the potential for detecting tumors at the earliest stages of 
development. Microwaves offer exceptionally high contrast compared to 
other imaging modalities^ such as X-ray mammography, which exploit 
intrinsic contrasts on the order of a few percent. Measurements suggest 
typical attenuation is less than 4 dB/cm up through 1 0 GHz, indicating 

that commercial microwave instrumentation with 100 dB of dynamic 

» 

range is capable of imaging through 25 cm of tissue. MIST uses 
microwave pulses that are on the order of 1 00 ps in duration, with peak 
powers on the order of a few milliwatts— approximately I/I OO^*^ of the 
power of a typical cellular phone. Assuming a pulse repetition frequency 
of 1 MHz and a maximum scan depth of 10 cm, an array of 100 antennas 
could be sequentiaiiy scanned in 0.1 seconds. 

[0025] The goal of conventional microwave tomography is the recovery 
of the dielectric-properties prof ile of an object from measurement of the 
transmission and scattering of hnicrowave energy through the object. In 
contrast, MIST in accordance with the invention need be carried out only 
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to identify the presence and location of strong scatterers in the breast by 
directly imaging backscattered signal power. Consequently, MIST avoids 
the need to solve a challenging, ill-conditioned nonlinear Inverse-scattering 
problem. Early active microwave backscatter techniques were 
unsuccessful because they used a single antenna location for transmitting 
and receiving and thus had no possibility of spatially focusing the 
backscattered signal. The use of an antenna array and short pulses 
enables MIST to focus in both space and time. Thus, MIST significantly 
enhances the response from malignant lesions while minimizing clutter 
signals, thereby overcoming challenges presented by breast heterogeneity 
and enabling the detection of lesions as small as 1-2 mm. Note that 
resolution is not determined by the wavelength of the microwave 
excitation. Rather, the spatial extent of the array aperture measured in 
wavelengths and the temporal duration of the pulse are the dominant 
factors in determining the resolution limit. 

[0026] Preliminary measurements suggest that the contrast between 
the dielectric properties of normal breast tissue and many benign lesions 
is negligible, in which case benign lesions would not act as strong 
microwave scatterers, allowing discrimination of benign and cancerous 
lesions. Furthermore, in contrast to conventional microwave tomography, 
MIST can also exploit morphology-dependent characteristics of lesions, 
such as spectral and polarization signatures, as vyell as the enhanced 
backscatter due to vascularization of malignant tumors, to further 
distinguish cancerous lesions from other scattering structures. In 
addition, change in lesion size is reflected in the backscattered spectral 
characteristics and signal-to-clutter ratio. 

[0027] An exemplary space-time microwave imaging system in 
accordance with the invention which provides transmission and reception 
with the same antenna is shown generally at 20 in Fig. 1 . The imaging 
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system 20 includes a microwave signal generator 21 whicli is supplied^ 
on a line 22, with clock pulses fronn a clock 23. The output of the signal 
generator 2 1 , which as described below may be short broadband pulses 
or a signal synthesized from multiple discrete frequencies, is provided on 
a line 25 to a power amplifier 26, the output of which is provided on a 
line 27 to a directional coupler 28. The output of the directional coupler 
28 is provided on a line 30 to a switching system 31 which selectively 
directs the power from the line 30 to lines 33 leading to each of the 
antennas 35 which are arranged in an array 36 of antennas (e.g., a 
rectangular or circular array). An array of antennas may be effectively 
provided by using one antenna 35 and moving it from position to position 
to collect data at each position, although the forming of a "virtual" array 
in this manner is not preferred- Further, the array may be formed to 
partially surround the object being imaged, for example, for use in breast 
imaging the array may formed to encircle the pendant breast. The 
antennas 35 and other microwave components should be wideband and 
preferably operate in the 1-10 GHz range. Examples of wideband antenna 
designs that may be utilized are the ''bowtie'' and Vivaldi type antennas 
and horn antennas designed for wideband operation. The switch 31 is 
formed to selectively provide a pulse of microwave power individually to 
the antennas 35 from the directional coupler 28 and to receive a signal 
from that antenna which is directed back through the switch 31 to the 
directional coupler 28. The directional coupler 28 sends the received 
signal on a line 38 to a low noise amplifier 40, the output of which is 
provided on a line 41 to a receiver 42. The receiver 42 also receives 
clock pulses on a line 43 from the clock 23. The clock pulses on the line 
43 allow the receiver 42 to time the onset of pulses of microwave power 
supplied from the signal generator 21 to allow correlation in time of the 
received signal with respect to the transmitted signal. Alternatively, the 
power output from the signal generator 21 may be provided through a 
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pawer splitter to the receiver 42 to allow time correfation. The signal 
generator 21, which may include a computer or digital processor, 
generates appropriately timed and shaped output pulses, discrete 
frequencies, etc., as required for the type of microwave transmission 
being utilized. The receiver 42 may be of conventional construction, 
providing detection of the received microwave signal and conversion of 
the detected signal to digitized data, e.g., with sampling of the received 
signal after each pulse to build up a digitized waveform, with the digitized 
data being provided to a digital signal processor of conventional design 
within the receiver 42 or to an appropriately programmed computer 44 
(e.g., a general purpose PC, a dedicated digital signal processor, etc.) all 
of which will be referred to herein generally as a "computer.'' It is 
understood that any type of computer that can be programmed to carry 
out the signal/data processing set forth herein may be utilized. The 
receiver 42 or the separate computer 44 processes the data to provide 
image data which may be displayed on a display device 45, such as a 
video display terminal, or which may be transmitted to a recording device 
46 such as a magnetic disk or CD ROM for long-term storageir or 
transmitted for printout, further data processing, etc. In accordance with 

■ 

the invention, space-time beamforming is carried out in a computer in the 
receiver 42 or a separate computer 44 on the data received from the 
antennas, as described further below. Further, signal processing is 
preferably employed to carry out a reflection artifact subtraction process 
(e.g., for the skin interface response or the antenna response) to reduce 
the effect of the artifact response on the received image data. Signal 
processing may also be carried out to compensate for frequency 
dependent scattering. As an example only of commercial instruments 
that may be utilized, the signal generator 21, amplifiers 26 and 40, 
directional coupler 28, receiver 42 and clock 23 may be implemented in 
an Agilent Vector Network Analyzer model 8720 ES, particularly for the 
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discrete frequency based approach, and the computer 44 may be 
connected to control the signal generator 21 and the switch 31 . 

# 

[0028] A space-time microwave irnaging system in accordance with the 
invention which may be utilized for simultaneous transmission from each 
antenna is shown generally at 50 in Fig. 2. The system 50 includes a 
signal generator 51 which receives a clock pulse on a line 52 from a clock 
53. The output of the signal generator 51 is provided on a line 54 to 
signal processing circuitry 55 which distributes the microwave (e.g., 
pulse) output on lines 57 to power amplifiers 58. Each of the power 
amplifiers 58 provides its output on a line 59 to a directional coupler 60, 
the output of which is provided on a line 61 to an individual antenna 63. 
The antennas 63 are arranged to form an array 64 of antennas, e.g., a 
rectangular array of antennas arranged in r-ows and columns. The signal 
processing circuitry 55 distributes the pulse of microwave on each of its 
output lines 57 with frequency dependent filtering to provide the desired 
microwave radiation from the antenna array 64, e.g., focussing of 
radiated power from the array 64 to selected points in the target object. 
The signals picked up by each antenna 63 are transmitted back on the 
line 61 to the directional coupler 60. The directional couplers provide the 
received signals on lines 65 to low noise amplifiers 66, the outputs of 
which are provided on lines 68 to a receiver 70. The receiver 70 also 
receives the clock pulses from the clock 53 on a line 71 to allow the 
receiver 70 to time the received signals with respect to the transmitted 
signals: The receiver 70 detects the microwave signal on a line 68 and 
converts the received signal to digital waveform data which is processed 
by a digital signal processor or a computer 72 in accordance with the 
invention. The image data from the computer 72 or digital signal 
processor may be displayed, e.g., on a video display terminal 73, or 
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provided to a storage device 74^e.g., CD ROM, magnetic disk, tape, etc. 
for long-term storajge, or transmitted for otiier purposes. 

[00291 Treatment for early-stage breast cancer typically involves a 
lumpectomy or partial mastectomy to remove the carcinoma and its 
margins, fallowed by radiation therapy to destroy any remaining cancer 
cells. For larger tumors, pre-operative chemotherapy may be used to 
shrink the tumor to conserve a larger portion of the breast. It is well 
known that the effect of radiation therapy and chemotherapy can be 
enhanced using microwave hyperthermia, that is, elevating the 
temperature of the cancerous tissue through microwave energy 
absorption. The persisting challenge in microwave hyperthermia, 
however/ is to preferentially heat the cancerous tissue without harming 
superficial and surrounding healthy breast tissues. Sophisticated adaptive 
focusing algorithms have been developed for use in phased-array 
hyperthermia treatment, but they require the use of invasive feedback 
probes located within the tumor. MIST technology offers a non-invasive 
approach for maximizing power deposition within the tumor and 
minimizing power deposition elsewhere. The microwave backscatter 
signals obtained during a low-power MIST scan of the breast inherently 
contain the information needed to tightly focus a transmitted high-power 
microwave pulse at the site of a tumor. In this manner, space-time 
microwave application utilizing the system of Fig. 2, configured to focus 
the microwave radiation at the position of the detected tumor, can 
provide hyperthermia treatment to destroy sifnall early-stage cancerous 
tumors without harming healthy tissue, thereby potentially eliminating the 
need for breast surgery and conventional radiation to the breast along 
with the accompanying side effects. 

[0030] With reference to Fig. 3, an antenna array device which may be 
utilized in the microwave imaging system of the invention is shown 
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generally at 80, having a face 81 over which are distributed multiple 
individual antennas 82 arranged in a two-dimensional array at known 
locations. The individual antenna elements 82 may have the "bow-tie" 
shape as shown or other shapes as desired. The array device 80 may be 
utilized as the antenna array 36 of Fig, 1 , with the antenna elements 82 
corresponding to the antennas 35, or as the antenna array 64 of Fig. 2, 
with the antenna elements 82 corresponding to the antennas 63. For 
purposes of illustration, the antenna array device 80 is also shown in Fig. 
3 placed adjacent to the breast 85 or other portion of the body to be 
imaged, preferably utilizing a matching element 86, such as a liquid filled 
bag, wliich conforms to the contour of the breast or other part of the 
body being imaged to minimize air gaps and unwanted reflections of 
microwave energy. While the invention is illustrated herein with regard to 
breast imaging, it is understood that the present invention may be utilized 
for imaging other parts of the body of an individual. 

[0031] To achieve the best resolution of the reconstructed image using 
the space-time focussing approach of the present invention, the radiated 
microwave pulse is preferably relatively short (e.g., about 100 ps), and 
thus has a wide band of frequency content, typically from 0 to 20 GHz 
and with significant energy in the frequency range of 1 GHz to 10 GHz, 
Thus, it is desirable to utilize antennas that are suitable for transmitting 
and receiving such short pulses with minimum distortion or elongation. It 
is desirable that the pulse radiating antenna have a constant sensitivity 
and a linear phase delay over the bandwidth of the incident 
electromagnetic pulse in the frequency domain. It is also desirable that 
the antenna design suppress both feed reflection and antenna ringing, and 
that the antenna have a smooth transition from the cabie impedance at 
the feed point to the impedance of the immersion medium at the radiating 
end of the antenna. The return loss, S11, should be low in magnitude as 
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less return loss means more power is transmitted to the antenna. Ideally, 
the return loss should be constant over the required bandwidth so that 
the spectrum of the transmitted power is flat and should have a linear 
phase delay across the frequency band so that the radiated waveform will 
not be dispersed. Other desirable properties include a weil-defined 
polarization, constant gain, and low side lobes in the radiation pattern. 
Resistively loaded cylindrical and conical dipole (monopole), and bow-tie 
antennas can be utilized for radiating temporally short, broad bandwidth 
pulses. Resistive loading can be utilized to reduce the unwanted 
reflections that occur along the antenna and the associated distortion of 
the radiated signal. Spiral antennas and log-periodic antennas have also 
been designed to achieve wide bandvyidth. Spectrum shaping and RF 
filtering may be needed to enhance the frequency performance of these 
antennas. Specialized antennas designed for pulse radiation may also be 

■* 

utilized;. An example of a suitable antenna that is designed for short pulse 
radiation is shown and described in U.S. Patent No. 6,348,898, issued 
February 19, 2002. 

[0032] As an example of the present invention, a MIST beamforming 
system was applied to simulated backscatter data generated from finite- 
difference time-domain (FDTD) computational electromagnetics 
simulations of microwave propagation in the breast. The anatomically 
realistic breast model was derived from a high-resolution 3-D breast MRI 
{magnetic resonance Imaging} obtained during routine patient care at the 
University of Wisconsin Hospital and Clinics. The face-down images of 
the pendant breast were digitally rotated, vertically compressed, and 
laterally expanded to create high-resolution images of the naturally 
flattened breast of a patient in a supine position. Then, each voxel was 
assigned the appropriate values of and a. The 2-D model is 

incorporated into FDTD simulations for a co-linear 1 7-element monopole 
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antenna array spanning 8 cm along the surface of the breast. Each 
antenna is excited with an ultrashort differentiated Gaussian pulse 
(temporal width of 1 1 0 ps, bandwidth of 9 GHz) and the backscattered 
response at the same antenna element is computed. This process is 
repeated for each element of the array, resulting in 17 received 
backscattered waveforms. The resulting FDTD-computed backscatter 
waveforms represent the scattering effects of heterogeneous normal 
breast tissue (clutter) and the malignant tumor (signal). 

[0033] The skin response subtraction process estimates the skin 
component of the signal at each antenna as a filtered combination of the 
signals at all other antennas. The filter weights are chosen to minimize 
the residual signal over that portion of the received data dominated by the 
reflection from an interface with the object being imaged such as the skin- 
breast interface. The results show that the skin response effect is 
removed at the expense of energy from the tumor bleeding throughout 
the image. This occurs because the skin response subtraction algorithm 
used somewhat distorts the response from the tumor. 

[0034] The beamformer algorithm utilized first time shifts the 1 7 
received signals to approximately align the returns from a hypothesized 
scatterer at a candidate location- The time-aligned signals are passed 
through a bank of finite-impulse response (FIR) filters (one in each 
antenna channel), summed, and time gated and the power calculated to 
produce the beamformer output signal. The filters are designed using a 
least squares tei^hnique to present unit gain to scattered signals 
originating from the candidate location. This technique is described in B. 
Van Veen, et aK, "Beamforming: A Versatile Approach to Spatial 
Filtering," IEEE ASSP Magazine, Vol. 5, April, 1988, pp. 4-24; B. Van 
Veen, "Minimum Variance Beamforming," \r\ Adaptive Radar Detection 
and Estimation t Ed. S. Haykin and A. Steinhardt, John Wiley and Sons; 
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New York, Chapter 4, March, 1992, pp. 161-236. Hence, the 
beamformer output power represents an estimate of the energy scattered 
by that location. The beamformer is scanned to different locations in the 
breast by changing the time shifts, filter weights, and time gating. The 
output power may then be plotted as a function of scan location/ Regions 
of large output power correspond to significant scatterers (e.g., malignant 

a 

lesions) . 

[0035] A simulation was carried out to determine the scanned MIST 
output power for a 2-mm-diameter malignant tumor located 3 cm deep. 
For this study, the average dielectric-properties contrast between 
malignant and normal breast tissue in the numerical breast phantom is 
approximately 5:1 . The heterogeneity of the normal breast tissue in the 
numerical breast phantom corresponds to variations in dielectric properties 
of ±10%, the upper bound on normal breast tissue variability that has 

been reported. The tumor was c/early detectable, as it stands out from 
the background clutter by 22 dB, MIST output power for two adjacent 2- 
mm-diameter tumors separated by 2 cm at a depth of 3 cm, showed two 
distinct scattering objects are clearly evident at the correct locations, 
demonstrating the potential resolving power of the present invention. A 
scenario under the worst-case assumption that the normal-tissue dielectric 
properties substantially exceed the published upper bound, thereby 
reducing the dielectric-properties contrast between malignant and normal 
tissue to less than 2:1, showed that even with significantly reduced 
contrast, the tumor was still easily detected, as the peak of the tumor 
response stands 1 1 dB above the largest background clutter. 

10036] The foregoing exemplary beamforming process incorporates 
frequency dependent propagation effects, but does not incorporate 
frequency dependent scattering effects. Scattering is frequency 
dependent due to dispersive dielectric properties and the presence of 

Mm 
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muttiple scattering surfaces. Frequency dependent scattering broadens 
the received pulse duration, reducing resolution, and slnifts the center of 
received energy in time, which causes scattered signal power to appear at 
an incorrect location. These errors may be compensated by processing 
the beamformer output signal from the filters prior to time gating using a 
parametric signal processing model for frequency dependent scattering 
effects- For example, autoregressive models may be used to describe the 
resonant behavior caused by finite tumor size. 

[0037] Removal of the response from the skin-breast interface is critical 
for lesion detection, as this response is orders of magnitude larger than 
the tumor response. This response may be removed at the expense of 
some distortion of the tumor response. The distortion is known, since it is 
a function of the weights used for skin response removal, allowing 
processing to be carried out for reducing or eliminating the tumor 
response distortion. 

^ - . 

[00381 The skin response removal algorithm estimates the skin 
response at each antenna. The skin response is a known function of the 
skin thickness and the dielectric properties of the skin and breast. This 
fact may be exploited in processes for estimating these properties from 
the skin response. The average breast dielectric properties may then be 
used as a calibration step to choose the best beamformer design for each 
patient, 

[0039] The methods described above assume only one antenna is 

■ ■ 

transmitting and receiving at any point in time. This process involves 
sequentially stepping through the array. If an antenna array with multiple 
receive channels is used as shown in Fig. 2, then a multitude of different 
trans mit-receive strategies are possfible. Beamforming and skin response 
removal algorithms may be utilized in which all antennas receive 

ft 
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simultaneously. Transmit strategies may also be utilized that focus the 
transmitted energy on a given region of the breast. The transmit and 
receive focus location is then scanned throughout the breast to form the 
image of scattered power. Such scanning may be utilized to improve 
resolution and robustness to artifacts, noise, and clutter. The signal 
parameters used to focus the transmission are the relative transmit time 
and signal amplitude in each antenna. After a lesion is located, the 
transmitted energy from the antennas may be focused on the lesion at a 
higher power level to heat and destroy the lesion. 

[0040] Methods may be employed for assessing changes in lesion size 
from images obtained at different points in time. Both the spatial extent 
of the scattering region as well as the total power returned may increase 
from one scan to the next if the tumor undergoes angiogenesis and 
growth. Tracking this growth would be useful in the diagnosis of . 
malignant lesions. Both the spatial extent of the scattering region and the 
total power returned may decrease if cancerous cells in the lesion are 
destroyed. Monitoring the decrease in lesion size would aid in assessing 
the effect of radiation therapy, chemotherapy, and/or thermotherapy . Use 
of absolute estimated tumor power is problematic due to expected 
variation from one measurement to the next. However, the peak tumor- 
to-clutter ratio should be robust to measurement variations and provide a 
reliable metric for assessing relative tumor size. Frequency dependent 
scattering effects will also vary with tumor size and provide another 
means for assessing changes over time. 

[0041] An exemplary MIST sensor in the imaging system of the 
invention may include a microwave vector reflectometer (the pulse 
generator 21, 51 and receiver 42, 70, and may include the associated 
amplifiers and directional couplers) and a low-reverberation ultrawideband 
transmitting/receiving antenna- A low-noise commercial vector network 
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analyzer (VNA) with a tEme-domain option may be used for the vector 
reflectometer. The dynamic range of a VNA of this type is sufficient to 
detect small malignant tumors up to depths of 5.0 cm in the breast. 

[0042] The MIST strategy for detection is to identify the presence and 
location of strong scatterers in the breast, rather than to attempt to 
reconstruct the dielectric-properties profile of the breast interior. As a 
resultjr the MIST approach overcomes the fundamental computational 
limitations and related vulnerabilities to noise of conventional narrowband 
microwave tomography. The use of spatial and temporal focusing in MIST 
significantly enhances the response from malignant lesions while 
minimizing clutter signals, thereby overcoming challenges presented by 
breast heterogeneity. Space-time focusing achieves super-resolutionr 
enabling the detection of extremely small (< 5 mm in diameter) malignant 
lesions with harmless low-power microwave signals. In contrast to earlier 
examples of breast imaging using ultrawideband microwave-radar 
techniques, MIST employs sophisticated and robust frequency-dependent 
processing of microwave backscatter signals to obtain superior sensitivity 
for discriminating against artifacts and noise* The innovative system 
configuration eliminates the need for breast compression and permits the 

r 

interior breast tissue to be imaged with the patient lying comfortably on 
her back. This uniquely enables MIST to detect tumors located near the 
chest wall or in the quadrant near the underarm where an estimated 50% 
of all breast tumors occur. 

[0043] Reflection artifact removal (such as skin response removal), 

T 

■ 

beam forming, and frequency-dependent scattering processes in 
accordance with the invention are discussed in further detail below. 
These processes may be carried out in a separate computer (e.g., the . 
computer 44 of Fig. 1 or 72 of Fig. 2), or in a digital signal processor of . 
the receiver (e.g., the receiver 42 of Fig, 1 or the receiver 70 of Fig. 2), 

21- 

■> 

BNSDOCID: <WO ^03003907A?_I_> 



wo 03/003907 PCT/US02/21397 

both of which will be referred to herein as a computer, that is 
programmed to carry out the processing on the digitized waveform signal 
data for each antenna that is provided by the receiver. 

[0044] The following describes the artifact removal and beamforming 
design method in mathematical expressions which are implemented in the 
computer and/or digital signal processors of the systems of Figs. 1 and 2. 
Lower and upper case boldface! Roman type is used to denote vector and 
matrix quantities, respectively. Superscript * represents the complex 
conjugate and superscripts T, H, and -1 represent the matrix transpose, 
complex conjugate transpose, and inverse, respectively. 

[00451 Reflection Artifact Subtraction 

[0046] A refiiection artifact removal process is preferably carried out on 
the data received from the antennas to remove large reflection artifacts, 
such as the energy reflected from the ends of the antenna and feed and 
from the sktn-breast interface. These reflections are typically orders of 
magnitude greater than the received backscatter signaL This reflection 
artifact removal or subtraction process will be described below for the 
example of removal of the skin-breast interface response. The skin 
response removal process forms an estimate of the response associated 
with the skin-breast interface and subtracts it from the recorded data. 
The response from the skin-breast interface is a function of the skin 
thickness, the dielectric properties of the skin, and the dielectric 
properties of the breast. Thus, the response from the skin-breast 
interface can be used to estimate these parameters. This is accomplished 
in general by expressing the response from the skin-breast interface as a 
parametric function of the unknown parameters and then choosing the 
unknown parameters to minimize the mean-squared error between the 
measured skin-breast response (the data) and the parametric function, 
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That is, wa choose the skin thickness and dielectric properties of the skin 
and breast so that the predicted response most closely approximates the 
actual measured response- 

[0047] To illustrate this, consider the simple case in which the dielectric 
properties are assumed to be frequency independent, the skin-breast 
interface is assumed planar, and the transmitted signal propagates as a 
plane wave. To further simplify the illustration/ we assume lossless 
propagation where the permittivity e^-, of the matching medium {e.g., the 

liquid matching medium 86 of Fig. 3) is known and the permittivity of the 
skin *he permittivity of the breast tissue s^g are unknown. Fig, 4 

illustrates the incident and reflected pulses at the skin/breast interface at 
four points in time, assuming an ultrashort microwave pulse. Fig. 5 
illustrates the observed temporal waveform at the antenna which includes 
the pulse reflected from the matching medium/skin interface (Ref #1) and 
the pulse reflected from the skin/breast tissue interface (Ref #2) (in 

J- 

practice, these reflected pulses will overlap). Let Vi^jclt) = v(t) represent the 

incident pulse. Then the first reflected pulse (Ref . #1) may be 
represented as 



t — 



[0048] and the second reflected pulse (Ref #2) may be represented as 



W -^12^23121^ 



t 



where 
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^12 * 



T = 



2^^ 



r2 



^r2 



and c is the free space propagation velocity. 

[0049] Hence, the observed waveform due to the incident pulse, the 
reflected pulse #1 and the reflected pulse #2 has the form 
r(t) -v{t) + a^v{t-t,) +a2v(t-t-,"t2) where v(t) is known, 

■ i' 

[0050] We may correlate v{t) with r{t) and estimate t,/ a, from the 
second peak (in time} and a^, from the third peak using standard time 
delay /amplitude estimation techniques. 

[0051] Given t^, a^, tgr we may solve for d-,, dg, ^>[s^ i y[s^ as 
follows: 

t 



where 
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[0052] There are many different techniques for estimating time delays 

1 

and amplitudes of a known waveform that are well known. The same 
general methods for determining the dielectric properties may be applied 
to more realistic models of the skin-breast interface. The following 
discusses the preferred solution of the skin response removal problem in 
further detail. 

[0053] Consider an array of N antennas and denote the received signal 
at the i^^ antenna as b^it}. Each received signal is converted to a sampled 

waveform, bjln], by an A/D converter in the receiver operating at a 

sampling frequency f^. The received signal contains contributions from 

the skin-breast interface, clutter due to heterogeneity in the breast, the 
backscatter from lesions, and noise. The response from the skin-breast 
interface is orders of magnitude larger than the response from all other 
contributions and thus must be removed prior to performing tumor 

detection- 

[0054] The skin artifacts in each of the N channels are similar but not 
identical due to local variations in skin thickness and breast heterogeneity. 
If the skin artifact for all channels were identical, one approach to remove 
it would be to subtract the average of the skin artifact across the N ^ 
channels from each channel. In order to compensate for channel to 
channel variation in the skin artifact> the skin artifact at each antenna may 
be estimated as a filtered combination of the signal at all other antennas, 
as shown in Fig. 6. The signals from each of the other antennas are 
provided to FIR filters 90, the outputs of which are summed at 91 and 
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subtracted at 92 from the signal from the particular antenna after a delay 
94, The filter weights of the FIR filters 90 are chosen to minimize the 
residua! signal mean-squared error over that portion of the received data 
dominated by the reflection from the skin-breast interface. Without loss 
of generality, suppose that the skin artifact is to be removed from the first 
antenna. Define the (2J -f 1 ) x 1 vector of time samples in the i**" antenna 
channel as 

h,[n]=[b,[n-Jl-'^b,[nl'-,b,[n^jf, 2<i<N (1) 

and let b2/^[w] = [b2 [«l -sb]^[nlf be the concatenation of data in channels 2 
through N. Similarly, let q,. be the (2J + 1) x 1 vector of FIR filter 

coefficients in the ^ channel and q = [qL"^q.^/]^ ^® concatenation of 

FIR (finite impulse response) filter coefficients from channels 2 through N, 
The optimal filter weight vector is chosen to satisfy 

q = argniiQ ^ h[«]-q'^*>2//[«] ^ (2) 

[0055] where n^ is the time that approximates when the skin artifact 

begins and m is the duration of the received signal that is dominated by 
the skin artifact. The solution to this minimization problem is given by 

q = R-'p (3) 



(4) 



[0056] The fact that there is a high degree of correlation among the 
skin artifacts in the N channels results in the sample covariance matrix R 
being ill-conditioned. If R is ill-conditioned, then the matrix inversion in 
equation (3) can result in a solution for q that has very large norm and 

■ 
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thus amplifies noise. In order to prevent this ^ we replace R with the low 
rank approximation 

(6) 

[0057] where A.•,^<\<p, are the p significant eigenvalues and U|,1<i<p, 

are the corresponding eigenvectors. The filter weight vector is 
determined by replacing R"' in equation (3) with 

K=tj-^^< (7) 

[0058] The skin artifact is then removed from the entire data record of 
the first channel to create artifact free data Xitnl given by 

^iW = ^>i[w]-q^b2j^[n] {8) 

[0059] This algorithm introduces a smaif level of distortion in the 
backscattered lesion signal because the backscattered lesion signals from 
the other N-1 channels are added back in to the first channel. This is 
explicitly shown by decomposing b,[n] and b^^In], into a skin artifact s,[nl 

and S2pj[n] and residuals d,[n] and d2N[n], respectively- The residual signals 

contain the backscattered response from the lesion. The values no and m 

are chosen so that q is determined from a portion of the data in which the 
residuals are negligible and, thus, 

[0060] However, decomposing bjn] and b^j^In] in equation (8) gives 

^d^[7t\-€i'd^^[n] (11) 
\ ■ -27' 
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[0061] Thus, the residual signal is distorted by q^dgwlnl. This term is 

generally small because q tends to "average" across channels and the 
lesion responses in d^^In] do not add in phase because they are not 

aligned in time. A simple method for reducing the distortion is to add a 
filtered version of the residual to obtain 

^1 [n] = ^, [»3 q^x^^ [n] (12) 
where 

^2N l^] = 1^2 [« - </! * • , ^2 ^ ' * ^ I" ~ -^1 * ^ ^ AT n 3 ) 

is the vector containing the data from the other N-1 channels after the 
skin artifact has been removed from each of them. This addition of a 
filtered form of the residual is illustrated in Fig. 1 1 which includes FIR 
filters 145 to provide filtered signals that are summed at 146 to produce 
a signal added at 148 to the corrected data signal x^[n\ to provide an 

improved corrected signal data x^[n]. - 

[0062] Fig. 8 are example waveforms showing the effect of the skin 
response subtraction process, with the solid lines indicating the original 
waveforms and the dashed lines indicating the waveforms after skin 
artifact removal. 

[0063] The artifact subtraction process can be applied only in the time 
domain. Thus, if frequency scanning is carried out using multiple discrete 

■ 

frequencies of the signals applied to the antennas, the received signal 
data must first be converted to the time domain (using an inverse FFT) 
prior to applying the artifact subtraction process. 

[0064] The artifact removal process requires that all of the artifacts 
occur at the same relative times in the different channels. If the antennas 
are located at varying distances from the skin, the skin response will 
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occur at different times. Thus, to apply the algorithm in general, the 
waveforms must first be time shifted so artifacts in all channels occur 
simultaneously. Aligning the artifacts in time is trivia) because by nature 
the artifact is huge and it is easy to see when it starts. 

[006f5] The antenna reflection response will not vary in time in the 
different channels (assuming nearly identical antennas), so time alignment 
is not needed for removing it. The algorithm can simultaneously remove 
antenna artifact and skin reflection artifact^ provided they are both time 
aligned in the waveforms. While this is true if the array is not the surface 
of the skin, it is not generally true if the distances to the skin differ for 
different antennas. In this case, one can apply the algorithm twice; first, 
to remove the antenna response, followed by time alignment of the 
residual skin response, and second remove the skin response. 

[0066] There is one limitation with applying it twice, and that has to do 
with the other requirement of the algorithm, which requires the artifact to 
be the only contribution to the signal over a time interval that spans at 
least part of the artifact duration. Hence, if the antennas are varying 
distances from the skin, but in some channels the skin response 
completely overlaps (in time) the antenna response, it may not perform 
adequately. 

[0067] Space-Time Beamforming 

[00681 The image of backscattered power as a function of a location r 
is obtained by scanning each location with a different space-time 
beamformer- The beamformer for scan location r forms a weighted 
combination of time-delayed versions of a signal as shown in Fig. 7, Each 
of the signals Xgln] from the antennas after subtraction of the skin 

response is passed through a time delay 97, a first time window 98, and 
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an= FIR filter 99. The outputs of the filters 99 are summed at 100- The 
summed output is passed through a second time window 1 02 to an 
energy calculation 103. The calculation of the energy in the signal over a 
consistent period of time provides a result which is proportional to power 
in the signal, and will be considered herein as the same as calculation of 
power. Preferably, the beamformer is designed to pass backscattered 

■ 

signals from the location r with unit gain, while attenuating signals; from 
other locations, 

10069] For design purposes, assume that the received signal on the i^^ 
channel is only comprised of the response due to a lesion at location r. 
Let this signal be denoted by Xj[n] having Fourier transform X,(ci>}. Note 

that the received signal is 

xX<&)^P{^yiif,a?) l<i<N (14) 

where P(o>} is the Fourier transform of the transmitted pulse p(t) and 
V|(r,a)) is the frequency response of the electromagnetic model 

representing frequency deperKJent propagation and scattering effects. 
The i**^ sampled waveform is then delayed by an integer number of 
samples njir) — ng-x,{r), resulting in the waveforms in each channel being 

approximately aligned in time. The average time Tj(r} denotes the 

roundtrip propagation delay for location r in the channel, computed by 
dividing the roundtrip path length by the average speed of propagation 
and rounding to the nearest sample, and n^ is the reference time to which 

all received signals are aligned. We choose n^ as the worst case delay 

over all channels and locations, that is, 

> round(pxax r (r )) (15) 
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[0070] The time aligned signals are windowed before the filtering stage, 
to remove interference and clutter prior to n^ that could contribute to the 

FIR filter outputs, using the window function 



g 



[n] 



1 if n>n„ 
0 otherwise 



[0071] The FIR filter in the i* channel has coefficients represented by 
the L X 1 vector Wj. The FIR filters equalize path length dependent 

dispersion and attenuation, interpolate any fractional time delays after 
time shifting, and bandpass filter the signal. The frequency response of 
each filter can be written as . 

FF;.(fi,)=gw„c--^"'^' =wfd(fl>) "(17) 

e-o 

where T„=1/i^ is the sampling interval and d{a>)=\l,e~^'^^' ,—,e-^'"^^^^^'Y . In 
order to pass signals from a location with unit gain and a linear phase 
shift, we require 

/=i M (18) 

where Fi(r^j,£») is the frequency response of the electromagnetic model 

after removing the average propagation and r^(L-l)/2 represents the 

average time delay introduced by the FIR filter. Thus, the design 
constraints on Wj are 

Y/,{^^,w)s^]a{€^^^ (19) 

[0072] If these constraints are satisfied, then the summed output of the 
FIR filter bank, z[n3, has the Fourier transform 
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Z(^xj) = P(^^?y^"^^^^^^^>'^^^^^ (20) 
[0073] This signal is windowed with h[rQ,n] to eliminate additional 

« 

clutter, and the power in the windowed signal is obtained by taking the 
sum of the square of each sample. The following sections discuss how 
the filter weights are designed and how the window h[rQ,n] is obtained. 

[00741 FIR Filter Design 

[0075] Let the NL x 1 filtering vector be w -[wf ,*--,w]^]^ * We may 
rewrite equation (19) as 

w^d(ro,£^?)«e-^"^-^^-^>^^ (21) 

where the NL x 1 array response vector corresponding to coordinate is 
d{r,,a))^d{o>)®[v,{r^^^^^ (22) 

[0076] Here the symbol ® denotes the kronecker product and 
superscript H denotes complex conjugate transpose. The filters are 
designed using a least squares technique to approximate equation (21) 
across the band of frequencies [o),r<oJ while attenuating signals from other 

directions. This is accomplished by seeking to approximate equation (21) 
on a dense grid of M distinct frequencies- To ensure that w is real- 
valued, vye use positive and negative frequency pairs. Define the matrix 
of array response vectors for M frequencies as 

A-[d(ro,a?,V -.^i-o.^^ii/)] (23) 
So, equation (21) is expressed as 

w^A«^f^ (24) 
where 
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Thus, the least squares design problem can be written as 



mm 

w 



. " 2 



(26) 



The minimum-norm solution to this problem is 



w = (AA^)"'Af^ (27) 

[00771 The solution may have a very large norm if A is ill-conditioned- 
A large norm can cause the gain at locations other than r^ to become 

large and also amplify noise. In order to control these effects, we choose 
w to solve the penalized least squares problem 



mm 



w 



' (28) 



2 



where X is a constant chosen to trade the norm of w against the 
approximation error. The solution is 

w = (AA^+;iI^^^J"'Af^ (29) 
[0078] Window Design 

[0079] The purpose of the second window 1 02 having the window 
function h[r^,, n] is to reduce residual interference and clutter in the signal 

zInL In the design process, the 2[n] is obtained by assuming the received 
signals are due to a point source at and by passing these signals 

through the designed space-time beamformer. The leading edge of the 
window is placed at the start of the beamformer output due to a point 
source. The response from finite size lesions have longer duratiorv than 
that of the point source due to frequency dependent scattering effects. 
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Hence, in order to capture the majority of the backscattered energy from 
the lesions, the window length should be chosen based on the 
approximate duration of actual expected responses from finite size 
lesions, if the window is too long, then the baclcground clutter level 
increases and resolution decreases. As an example, the length of the 
window may typically be selected for tumor sizes having a diameter in the 
range of 2 mm to 8 mm. 

[0080] Frequency Dependent Scattering 

[0081] A lesion may Introduce frequency-dependent scattering due to 
frequency dependence in its dielectric properties (and those of the 
surrounding medium}. There is also the effect of energy resonating in the 
"cavity" represented by the lesion (reverberant scattering within the 
lesion). This effect will depend on the size of the lesion and its 
morphology. 

[0082] We expect frequency dependence in the scattered signal to be 
of use in monitoring lesion growth as well as diagnosis (e.g., 
differentiating malignant from benign based on morphology). . 

[0083] The model is v 



[0084] This assumes a "point scatterer" or frequency-independent 
scatterer since P{oj) is the Fourier transform of the transmitted pulse anff 
Vi(r, o) Is the propagation path frequency response. 

[00851 A first-order frequency-dependent scattering model is 



xXg>) P{a>)V,{Y,a)) \<i<N 



(30) 



xM^P{^)A<^)Vi^o^o>) \<i<N 



(31) 
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[0086] This model assumes the lesion modifies the spectrum of the 
incomrng pulse with A{co). This is a first-order model because it ignores 
possible dependence on i. 

[0087] The beamformer output for location (prior to time gating) has; 
the frequency response 

Z(ro,£i?) = P(^i?)4ro,^)r(ro,a>) (32) 
where 

T{r,,o^) = f;^VX^,,a>)wfd{^^^^ (33) 

is the frequency response due to propagation and beamforming, and is 
known given r^r Wj, n^, and the average dielectric properties. P{ce)) is also 

assumed known (or can be measured). The problenri is to identify A(ro,o)). 

This is a classic system identification problem 

[0088] We have A(r^,w)-= ^^^^f^ 

[0089] The division in the frequency domain shown here is for 
conceptual purposes. In practicei. this approach is sensitive to noise and 
small values of Pico) T(ro, o) and other system identification methods nfiay 

be used. 

[0090] It is in general preferable to use a parametric model for A{ro, O)) 

to minimize modeling of extraneous clutter at the output of the 
beamformer. Standard models that may be used to this end include 
autoregressive (all pole)^ moving^average {all zero), and autoregressive- 
moving average (pole-zero). Methods for determining model parameters 
from system input and output data are well known. 
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C0091] To illustrate, consider using an all-pole (AR) model. This is wei! 
suited to describing resonance effects. We shall formulate this problem in 
the time domain. Let pj{n) be the beamformer output calculated assuming 

a point scatterer, and zlr^^n) be the actual output. We drop from now 

on . 



[0092] Here a: represents the AR coefficients. The model states 



aoZ(n) + a^z{n-1 ) +a2z(n-2) + ... 4-a 2{n-p} = PT-(n) 



(34) 



[0093] The presence of noise and modeling error will cause this to be 
an approximation. One method for choosing the aj's is to minimize the 

approximation error. 

[0094] Given z(n), pj{n)> n = 0,1 ,-.Q (Q>P) we may write this as 
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(35) 



[00951 One means for identifying a is to solve this problem in the least- 
squares sense. Many other methods can also be formulated, as is well 
known. 

[0096] The frequency response is expressed in terms of the a/s as 



(36) 



[0097] Given a model for A(o)}, one can potentially discriminate 
between classes of lesions or lesion size, it is also helpful to "equalize" 
the effects of the frequency-dependent scattering prior to time gating the 
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beamf ormer output and forming an image. This equalization can reduce 
the sensitivity of the time-gate duration to tumor size. The duration of 
the backscattered signal varies due to frequency dependent scattering and 
thus if these effects are not equalized, the time-gate duration must be 
. chosen long enough to capture sufficient energy from lesions with long 
duration backscatter. This reduces resolution for lesions with short 
duration backscatter. Thus, a preferred method for accommodating 
frequency-dependent scattering is shown in Fig. 9, As iilustratad therein, 
the output of the beamf ormer 110 is provided to an equalizer 111 and to 
a system identifier 112 (which also provides control signals to the 
equalizer 111). The output of the equalizer 111 is provided through a 
time gate 1 1 4 to a power calculator 115 which provides its results to an 
image generator 117, Information from the image and from system 
identification is provided to a lesion discrimination/properties estimator 
11 8. 

[0098] Various techniques for equalizing signals are well known in the 
communications and signal processing literature and may be utilized- 
Note that the equalizer should not remove the average "gain" of the 
frequency dependent scattering sincp this factor is needed to get a 
correct relative estimate of backscattered power. 

[0099] Beamformer Design in the Frequency Domain 

[0100] A frequency-domain design approach may also be utilized for 
microwave imaging vf a space-time beamforming. The frequency domain 
design MIST beamformer spatially focuses the backscatter waveforms, 
compensates for frequency-dependent propagation effects, and 
discriminates against artifacts and noise. This is accomplished by time- 
aligning the backscatter, followed by frequency dependent filtering in 
each antenna channel. The filter frequency responses preferably solve a 
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penalized least squares problem designed to pass signals from a candidate 
scattering location in the breast with unit gain, subject to soft constraints 
on the norms of the filters for each channel. The MIST beamforming can 
produce a 3-D image of the relative backscatter energy, where each pixel 
of the image represents the energy of the backscatter originating from the 
corresponding location in the breast. Image locations with high energy 
implicate the presence of malignant tissue. The frequency-domain MIST 
design can provide a less computationally Intense solution. Given M 
channels and L filter coefficients per channel, time-domain MIST 
beamformer design requires inversion of an ML x ML matrix for each pixel 
in the design stage. The frequency-domain design circumvents matrix 
inversions altogether at the cost of selecting multiple penalty weights 
rather than a single one. 

[0101] The MIST beamformer for either the time- or frequency-domain 
designed filters may be implemented in the frequency domain using a 
series of single frequency measurements to synthesize a broadband input. 
Assuming M channels and A/ DFT {discrete Fourier transform) frequencies, 
the imaging procedure involves an /V-point FFT {fast Fourier transform) on 
the data in each of the M channels, one IFFT (inverse fast Fourier 
transform) plus the matrix multiplication of two N x M matrices. Thus, 
the number of computations for each Image pixel is on the order of 

[0102] For simplicity we consider the monostatic case which is a 
single-transmit, single-receive configuration. Extensions to the single- 
transmit, multiple-receive configuration are straightforward. In the 
monostatic case each antenna, in turn, transmits a pulse into the breast 
and collects the resulting backscatter before the next antenna repeats the 
sequence. The collected signals are converted to sampled waveforms 
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containing contributions from the skin-breast interface, clutter due to 
heterogeneity in the breast, backscatter from possible lesions, and noise. 
Since the response from the skin-breast interface is much larger than the 
response from all other contributions, it must be removed prior to 
performing tumor detection . We accomplish this with the skin-breast 
artifact removal algorithm presented above which estimates the skin- 
breast artifact in each channel as a filtered combination of the signals in 
all other channels. The filter weights are chosen to minimize the residual 
signal mean-squared error over that portion of the received data 
dominated by the reflection from the skin-breast interface. This algorithm 
is shown to effectively eliminate the skin-breast artifact while introducing 
only a very small level of distortion in the signal faackscattered from the 
lesion. 

[0103] The image of backscattered energy as a function of scan 
location r is obtained by applying a space-time beamformer designed for 
each scan location to the backscattered signals. Fig. 10 illustrates the 
post-processing performed on the signals for a specific scan location, r^, 

following the skin-breast artifact removal. The signals X,-[nl are time 

aligned with an appropriate time delay 1 30, windowed at 132, 
transformed to the frequency domain at 134 and pointwise multiplied by 
the beamformer coefficients at 135. The sum of these weighted signals 
at 137 forms the beamformer output An Inverse DFT 138 

transforms the beamformer output back to the time domain where a time- 
gating window is applied at 139* Finally, the time-gated signal energy is 
calculated at 1 40, resulting in the backscatter strength assigned to 
location in the image. 

[0104] For design purposes we assume that the received signal in the 
z^' channel contains only the backscatter due to a lesion at location . 
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Let this received signal after skiivbreast artifact removal be denoted by 
X'[n] and its DFT denoted 

X,[£]^l{a},)Si,{r^,a},\ l^i^M,l<£^N (37) 

where a)g is the frequency corresponding to the £'^' DFT index, l{ci}^) is 

the DFT of the transmitted pulse and S^^ir^^co^) is an analytical model of 

the monostatic frequency response associated with propagation through 
breast tissue from the antenna to the scatterer located at and back. 

Note that this is the same relation as Eqn. (14) except for the use of 
discrete frequencies. l{co^) corresponds to P{o)) in Eqn. (14) and 

*S'i,(?o>^^) corresponds to I^(r,c2?) in Eqn. (14). This notation is used below 

to illustrate the beamformer design in the frequency domain. A series of 
narrowband measurements rhay be used in place of a single 
ultrawfdeband measurement for each antenna location, in which case 
represents the amplitude and phase of the transmitted 

sinusoidal signal with carrier frequency 

[0105] As in the time domain design, we delay the signal jcJ/z] at 130 
by an integer number of samples n^{r^)^n„ ""^/O^o) so that the waveforms 
in each channel are approximately aligned in time. Here is the 
reference sample interval to which all received signals are aligned and 
r,(ro) denotes the round-trip propagation delay for location in the 

channel in units of sample intervals* The delay is computed by dividing 
the round-trip path length by the average speed of propagation and 
rounding to the nearest sample. We choose as the worst case delay 

over all channels and locations, that is, 

»,=maxi:,(rj (38) 
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[01061 The time-aligned signals are windowed at 132 before the 
filtering stage to remove interference and clutter that is present prior to 
time n. using the window function 



S 



1 n>n„ 

(39) 

0 ofhenvise 



At this stage the signals are transformed to the frequency domain at 134 
and passed through the beamformer at 1 35. The purpose of 
beamforming is to equalize the transmitted pulse, equalize path-length 
dependent dispersion and attenuation, interpolate any fractional tirne 
delays, and perform spatial discrimination against clutter. The 
beamformer weight in the z'^ channel at DFT frequency index t is 
denoted as Wi\i\. We design the beamformer using the minimum number 

of coefficients or equivalently the minimum number of design frequencies. 
The minimum number of frequency-domain MIST beamformer 
coefficients, L, may be determined as follows. The maximum distance of 
any point in the scan region from the antennas determines the time extent 
of interest for the received signal, and consequently determines the 
minimum DFT length, M The DFT length, N, and the sampling frequency, 
, determine the spacing between the DFT frequencies as Aa} = a}-/N. 

If the bandwidth of interest is B , then the minimum number of 
frequencies that must be considered in the beamformer design is given by 

BN 



(40) 



We assume that the band of interest corresponds to the set of £ 
frequencies, {o)^ :£q ^1^£^+L-1}. 

In order for the beamformer to compensate for the transmitted pulse and 
propagation to and from location r^,, we require the output to have unit 
gain and linear phase as shown by 
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101071 Here S^ij^Q^co^) denotes the frequency response due to 

propagation after removing the linear phase shift associated with round- 
trip propagation delay, r^ir^), the value t^={N-\)12 represents the 

average time delay introduced by the beamformer, and is the sampling 

interval. Combining the phase factors associated with the propagation 
and the time alignment according to ^ r.{r(j)+Mf(r^), we obtain the 

design constraints on J?^-[/] as 

[01081 Stack the Wi{i\ and 5a(r„,<»^) to form the Mxl vectors 

■ • '"'^ 

S(ro,a>f) = /(iyj[^a(ro,iOf)-S22(ro»<»J"^^ (43) 

^{^>\^Mw^{i\^'^W^{&J (44) 

and rewrite equation (42) in compact form as 

W^[^p(ro.G>,) = e--''"''«^S i^<,i<£^+L-\ (45) 

[0109] If these constraints are satisfied, then the beamformer output at 
frequency is given by 

Z(aj,) = e"^"'t'^^"'-^''' i^<i<Jt^^L-\ (46) 

[0110] The beamformer output is converted back to a time-domain 
signal^ z\fj\, using an inverse DFT at 138 and a window, ^[I'o^^w]/ is applied 

at 139 to eliminate additional clutter. The output energy at this scan 
location, p^f^)* is obtained at 140 from the sum of the squares of the 

windowed signal 
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Pi^o) = T,Hn]hKn]f _ . (47) 



n 



[0111] The reconstructed image of microwave scattering strength is 
obtained by scanning Xq throughout the reconstruction region and plotting 

beamformer output energy as a function of location. 

[011 2] The following subsections describe details specific to the design 
of the filter weights w[£] and the design of the window A[ro,K]- 

[01131 Beamformer Design 

[01 14] Designing the beamformer in the frequency domain has the 
advantage of yielding an exact solution to Eqn. {45) which can be 
expressed in closed-form: 

[01 15] This solution is especially appealing in comparison to the time- 
domain approach since it requires no matrix inversion. However, the 
exact solution is not robust because the magnitude of the beamformer 

weights can become very large when S^(ro,fi?^)s(ro,ii?^) is small. This is 

particularly problematic at higher frequencies and deeper scan locations 
where attenuation results in smafi values for S^Itq^co^). The robusitness of 

a beamformer to errors between actual and assumed propagation models 
and to background noise is proportional to the norm of the weight vector 
or noise gain. Using Eqn, (48), we have the noise gain 

gJ^]=W^MwM= ^^ = ^ _ (49) 
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[0116] where the subscript Is refers to the least squares (exact) 
solution for W[i]. Thus> the noise gain of the exact solution can become 

arbitrarily large as S^.(ro,fi?^)j decreases- 

[0117] In order to control the noise gain and obtain a robust 
beamformer, we consider the penalized least squares problem 



2 



(50) 



[0118] where the first term is the approximation error, the second term 
is the penalty function, and the penalty weight, is selected to trade 

the norm of against the approximation error. Defining a[^] as an 

■ 

MxM diagonal matrix with elements the solution to 

Eqn. (50) may be written as 

l + S^(ro,a?,)A-^[-^p(ro,a?,) 
[0119] Note that as /L,-[^] decreases the approximation error decreases 
but the norm of W^li] increases. In order to compromise between these 

effects, we select the penalty weights >l|[^] = [s,y(ro,<2)f)| . With this choice 

the beamf ormer coefficients simplify to 

[0120] This solution for the beamformer weights does not require 
matrix inversion. 

[0121] We may compare the approximation error and noise gain of this 
penalized least squares beamformer of Eqn. (52) to those of the 
unconstrained beamformer in Eqn. (48). The approximation error of Eqn. 
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(48) is zero and its noise gain is given by Eqn. (49), For the beamformer 
of Eqn. {52}^ the approximation error is given by 



and the noise gain is 



where the subscript p/s refers to the penalized least squares solution for 
Comparison of Eqn, (49) and Eqn. (54) reveals that 

<jp& M ^ [^] - Furthermore, Eqn, (49) indicates that the noise gain of 

the penalized least squares beamformer cannot exceed M. 
[0122] Window Design 

[0123] If the beamformer satisfies Eqn. (45) and the lesion is a point 
scatterer, then the output z[nj is a time-shifted, attenuated and sampled 
version of a bandlimited impulse. If the values of z{ca^) outside the band 

of interest are set to zero, then z[nl has a sine function envelope and the 
majority of backscattered energy is contained within the mainlobe. Since 
both the time shift and bandwidth are known, the location and width of 
the mainlobe are also known. If the mainlobe occupies time points 

through mz[n], then a natural choice for the window is 

0 otlierwise 

[0124] This choice reduces clutter effects by ensuring that the output 
energy of Eqn, (47) is calculated using only samples of z[n] containing 
backscattering lesion energy. 
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[0125] In practice, scattering from the tumor is frequency-dependent, 
so the beamf ormer output corresponds to the impulse response of the 
tumor on the band of interest- These dispersive effects increase the 
duration of the beamf ormer output and complicate window selection. Our 
preliminary investigations suggest that the extent of the increase in 
duration is directly proportional to the tumor size. Since we are interested 
in detecting very small lesions, we have chosen to design hlr^rn] 

assuming a point scatterer model. This gives the largest possible signal- 
to-cfutter ratio (S/C) for small tumors. The S/C for large tumors is 
reduced by this choice; however, the backscattered signal from larger 
tumors is much stronger so a compromised S/C is relatively 
inconsequential for tumor detection. 

[0126] Example of a 2-D Spacer-Time Beamformer Design 

[0127] In order to tllustrpte the MIST beamf orming algorithm presented 
above, we discuss the frequency-domain design of a 2-D space-time 
beamformer. For ease of presentation, we show only the 2-D case, but 
these techniques are directly applicable in three-dimensions (3-D). 

[0128] Frequency-dependent propagation effects are incorporated into 
the design of the space-time beamformer via a monostatic transfer 
function, 5'^(r,ftj), which relates the received signal at the antenna at 

location to the transmitted signal at the antenna due to a scatterer 

located at r. In our 2-D example, we assume that each antenna is an _ 
infinite line source of electric current, I{cd), located atr{ in a uniform 

medium of normal breast tissue. The scatterer is treated as a conducting 
circular cylinder of infinitesimal radius and infinite length. For this 
scenario, we use the following transfer function: 
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1 



(56) 



where aipo) is the frequency dependent attenuation constant and fi(ai) is 
the frequency dependent phase constant. 

[0129] The 2-D beamfprmer is designed for a rigid 1-D conformal 
antenna array. The array contains 1 7 elements spanning 8 cm horizontally 
along the surface of the breast at known locations. The 2-D place that 
the beamformer is designed to scan spans 10 cm and is 4 cm deep. The 
transmitted UWB pulse is a differentiated Gaussian with a full width at 
ha!f maxihnum equal to 110 ps. Its spectrum has a peak near 6 GHz and 
significant energy between 1 and 1 1 GHz. 

[0130] We design the frequency-domain beamformer over the band 0.5 
to 12 GHz and apply a raised cosine taper over the intervals 0.5 to 1 GHz 
and 1 1 to 1 2 GHz to smooth the transition between passband and 
stopbands. The sampling frequency, /g, is 50 GHz and the maximum 

round-trip distance between any antenna and the deepest scan location is 
24 cm which corresponds to a time interval of /v^^ 125 sampling intervals. 

Thus A/= 1 25 is the data record length and L = 28 frequencies are 
employed to span the band of interest. The design location r is scanned 
over the breast region using a grid resolution of 1 mm. The post- 
beamformer window described by (19) is three sampling intervals in 
length, spanning 60 ps. 

[0131] After applying the skin-breast artifact removal algorithm to the 
simulated data, the space-time beamformer is applied to each scan 
location in the breast. The signal energy is approximated by steering the 
beamformer to the center location of a 2-mm-diameter tumor, applying 
windows of varying length to the beamformer output, and calculating the 
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output energy- Similarly, the clutter energy is computed by applying the 
same process to backscatter waveforms obtained from a tumor-free 
model, repeating for each scan location within the central 6 cm x 2 cm 
region of the breast, and averaging the results. We choose a window 
length of three samples (60 ns) in order to balance the goal of maximizing 
S/C with the goal of maintaining robustness to dispersion. A larger 
window would capture more signal energy at this scan location, but that 
same window would capture more clutter energy when the beamformer is 
steered elsewhere in the breast. The goal of preserving signal energy 

■ 

while discriminating against clutter suggests this strategy of optimizing 
S/C rather than signal energy alone. 

[0132J Transmitted Microwave Focussing 

[0133] As indicated above, the transmitted microwave energy from the 
antennas may be focussed at a chosen location for purposes such as 
increasing the signal-to-clutter ratio and {at sufficiently high energies) 
hyperthermia treatment of a tumor that has been located as described 
above. The appropriate timing and shaping of the pulses from the various 
antennas may be carried out in the signal processing circuitry 55 of a 
system as shown in Fig. 2, The configuration of the signal processing 
circuitry 55 for this purpose (which may be implemented in a digital signal 
processor) is illustrated in Fig. 1 2. The microwave pulses from the signal 
generator 51 are provided on the lines 54 to a time delay unit 1 52 for 
each of the N antennas 63 in the array 64. The outputs of the delay units 
1 52 are provided to FIR fitters 1 55, and the filtered signals are provided 
on the lines 57 to the power amplifiers 58. The following describes the 
design of the delays 152 and filters 155, 

[0134] Let the propagation model from the i*^ antenna 63 to the 
hypothesized tumor location (r^) be {t^yO}) and (r^) be the one-way 
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propagation delay. As before, we choose «_ as the worst case delay 



over ail channels, that is, n„ = 



max 



Tiir^) . Assume r J. (To) is rounded to 



the nearest integer of samples. Let {r^,cD) = V^ (r^^ayje^'^^'^^^'K Writing 
W,{oi}) = wJd{ct}) as in Eqn- (17), the design equations for the FIR filters 155 

are ^f^. (ro,<2>Vfd(^y)» e"^**"^^^^^ is analogous to Eqn. (19). Solve 



1=1 



this as in Eqn. (28) (penalized least squares sense), with the obvious 
changes to Eqn- (22). 

[0135] We may pick the tumor location from the peak of the irhage. in 

the presence of errors between the assumed propagation model (fo^^) 

and the true one, it is probably best to choose the peak from an image 
generated by transmitting from the center antenna and receiving at all 
antennas, since then the difference between propagation paths to each 
antenna is the one-way trip (the trip from antenna to tumor is the same in 
ali receive channels). In this way, errors in the propagation mode/ that 
lead to localization errors should be approximately compensated. 



10136] Images generated by MIST beamforming do not uniquely 
determine the exact tumor configuration in the breast. Deducing the 
tumor configuration from the energy plots does not have a unique solution 
because different tumor configurations couid lead to qualitatively similar 
energy plots. If, however, on the assumption that at most a single turnor 
is present then the image does provide insight into tumor properties, such 



[01371 The foregoing examples of beamformers were designed for the 
actual average dielectric properties of the normal breast tissue present in 
the model. In practice, these properties will vary from patient to patient 
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within a certain margin, and exact normal breast tissue properties will be 
unknown. The frequency-domain design is nonetheless robust with 
respect to potential mismatch between the assumed and the actual 
dielectric properties averages. 

% 

[0138] It is understood that the invention is not limited to the 
embodiments set forth herein for purposes of illustrating the invention, 
but embraces all such forms thereof as come within the scope of the 
following claims. 
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CLAIMS 



WHAT IS CLAIMED IS: 

1 1 . A space-time microwave imaging system for imaging 

2 an object comprising: 

3 (a) an array of antennas for radiating and receiving 

4 microwaves; 

5 (b) a microwave pulse source connected to tiie 

6 array of antennas to provide pulse signals of a selected width and 

7 repetition rate to the antennas; 

8 (c) a receiver connected to the antennas to detect 

9 the microwave signals received by the antennas and provide signal data 

10 corresponding thereto; and 

11 (d) a computer connected to receive the signal 

12 data, the computer programmed to estimate an artifact reflection 

13 component of the signal at each antenna as a filtered combination of the 

14 signals at ail other antennas and to subtract the estimated artifact 

15 reflection component from the signal data to provide corrected signal 

16 data. With weights of the filters chosen to minimize the residual signal 

17 over that portion of the received data dominated by the reflection, and the 

18 computer further programmed to process the corrected signal data in a 

19 beamformer process to time shift the corrected signal data to 

20 approximately align the returns from a scatterer at a candidate location, to 

21 pass the time aligned signals through a bank of filters with a filter for 

22 each antenna, to sum the outputs of the filters to form a summed signal, 

23 and to calculate the power in the summed signal to produce a 

24 beamformer output signal, wherein the beamformer filters are designed to 

25 present maximum gain to scattered signals originating from the candidate 

26 location, and to scan the beamformer process to a plurality of different 
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27 candidate locations in the object to be imaged by changing the time shifts 

28 and filter weights to generate multi-dimensional output power data. 

1 2- The imaging system of Claim 1 including an output 

2 device connected to the computer to display the raulti-dimenslonal output 

3 power data as a function of scanned locations. 

■ 

1 3. The imaging system of Claim 1 wherein the computer 

2 is further programmed to time gate the summed signal to form a time 

3 gated summed signal and to calculate the power in the time gated 

4 summed signal, and to scan the beamformer process to a plurality of 

5 different candidate locations in the object to be imaged by changing the 

6 time shiftSr filter weights and time gates to generate multi-dimensional 

7 output power data. 

1 4. The imaging system of Claim 3 wherein the computer 

2 is further programmed to process the beamformer output signal from the 

3 filters prior to time gating using a parametric signal processing model to 

4 compensate frequency dependent scattering effects. 

1 5, The imaging system of Claim 1 wherein the 

2 microwave pulse source provides pulses having pulse widths on the order 

3 of 100 picoseconds or less in duration* 

1 6. The imaging system of Claim 1 wherein the 

2 microwave pulse source is connected to the antennas to provide pulse 

3 signals to one antenna at a time. 

t 7. The imaging system of Claim 1 wherein the 

2 microwave pulse source is connected to the antennas to provide pulse 

3 signals to all of the antennas simultaneously. 
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1 8- The imaging system of Claim 1 wherein the computer 

2 is programmed in the beamformer process to apply a selected window to 

3 the time aligned signals before passing the time aligned signals through 

4 the bank of filters, and further to apply a selected window to the summed 

5 signal before the power in the summed signal is calculated. 

1 9. The imaging system of Claim 1 wherein the computer 

2 is programmed to calculate the power at each candidate location as the 

3 Slim of the squares of each sample of the summed signal. 

1 1 0. The imaging system of Claim 1 wherein the 

2 beamformer filters are FIR filters designed to satisfy a penalized least 

3 squares condition to present unit gain to scattered signals originating from 

4 a candidate location. 

1 11. The imaging system of Claim 1 further including signal 

2 processing circuitry that receives the pulses from the source and passes 
a the pulses through a delay and a filter for each antenna before providing 

4 the delayed and filtered pulses to the antennas, the delays and filters 

5 selected to focus the radiated microwave energy from the array of 

6 antennas at a selected candidate location in the object. 

1 12. A method of carrying out space-time microwave 

2 imaging of an individual comprising: 

3 (a) transmitting microwave signals from a plurality of 

4 antenna locations into an Individual to be examined; 

5 (b) receiving backscattered microwave signals at a 

6 plurality of antenna locations to provide received signals from the plurality 

7 of antenna locations; 

8 (c) processing the received signals in a computer in a 

9 beamformer process by time shifting the received signal data to 

10 approximately align the returns from a scatterer at a candidate location, 
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11 passing the time aligned signals through a banic of filters with a filter for 

12 each antenna location, summing the outputs of the filters to provide a 

13 summed signal, and time-gating the summed signal to provide a time- 

14 gated signal, and calculating the power in the time-gated signal to 

15 produce a beamformer output signal, wherein the filters are designed to 

16 present maximum gain to scattered signals originating from the candidate 

17 location; and 

18 (d) then scanning the beamformer process to a plurality of 

19 different candidate locations in the individual and repeating steps (a), (b) 

20 and (c) at each candidate location by changing the time shifts, filter 

21 weights and time gates for each candidate location to generate multi- 

22 dimensional output power data. 

1 13- The method of Claim 12 wherein the step of 

2 transmitting microwave signals comprises transmitting microwave pulses 

3 having pulse widths on the order of 100 picoseconds or less in duration. 

1 1 4. The method of Claim 1 2 wherein the step of 

2 transmitting microwave signals comprises transmitting microwave pulses 

3 having frequency content at 1 0 GHz or higher. 

1 15. The method of Claim 12 wherein the step of filtering is 

2 carried out in the time domain . 

1 16- The method of Clajm 1 2 wherein the step of filtering is 

2 carried out in the frequency domain. 

1 17- The method of Claim 1 6 wherein the step of filtering 

2 in the frequency domain is carried out for each signal from an antenna 

3 location in a channel by perforrning a fast Fourier transform on the signal 

4 data to provide frequency domain data, applying weights to the frequency 
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5 domain data in each channei and^adding the weighted data; and then 

6 performing an inverse fast Fourier transform on the weighted data to 

7 provide a filtered time domain data signal. 

1 .18- The method of Claim 1 2 including compensating the 

2 received signal data for frequency dependent scattering before carrying 

3 out faeamforrner processing. 

1 19- The method of Claim 1 2 including transmitting the 

2 microwave signals from ah array of antennas so as to focus the 

3 microwave power on a candidate location. 

1 20. The method of Claim 12 further including, before the 

2 step of processing the received signals in a beamformer process, the step 

3 of: 

4 carrying out an artifact response subtraction process on the 

5 received signals in the computer by estimating the interface reflection ' 

6 component of the signal at each antenna location as a combination of the 

7 received signals at the other antenna locations passed through filters to 
Q provide corrected signal data, the filters having weights chosen to 

9 minimize the received signal over that portion of the received signal 

10 dominated by the artifact response, and providing the corrected signal 

11 data to the beamformer process. 

1 21, The method of Claim 1 2 wherein microwave signals 

2 are provided to one antenna at a time and backscattered microwave 

3 signals are received from one antenna at a time for each of the antenna 

4 locations. 



-55- 

BNSOOCID: <WO^_^_03003©07A2 I > 



wo 03/003907 PCT/US02/21397 



1 22. The method of Claim 12 wherein microwave signals 

z are transmitted from all of the antennas simultaneously and backscattered 

3 microwave signals are received from all of the antennas simultaneously. 

1 23- The method of Claim 12 including in the beamformer 

2 process the steps of applying a selected window to the time aligned 

3 signals before fiitering the signals and further applying a selected window 

4 to the summed signal before calculating the power in the summed signal. 

1 24. The method of Claim 12 wherein the power is 

2 calculated at each candidate location as the sum of the squares of each 

3 sample of the summed signal. 

1 25. The method of Claim 12 wherein filtering is carried out 

2 utilizing FIR filters designed to satisfy a penalized least squares condition 

3 to present unit gain to scattered signals originating from a candidate 

4 location. 

1 26. The method of Claim 12 wherein the step of 

2 transmitting the microwave signals is carried out simultaneously from aii 

3 of the antenna locations by passing microwave pulses for each antenna at 

4 an antenna location through a delay and a filter for each antenna, the 

5 delays and filters selected to focus the radiated microwave energy from 

6 the antennas at a selected candidate location in the object, 

1 27. A method of compensating for an artifact reflection in 

2 microwave imaging of an individual comprising: 

I- 

3 (a) transmitting microwave signals from a plurality of 

4 antenna locations into an individual to be examined; 



-56- 



03003907A2 I > 



wo 03/(103907 PCT/US02/21397 



5 (b) receiving backscattered microwave signals at a 

6 plurality of antenna locations to provide received signals from the plurality 

7 of antenna locations; 

8 (c) carrying out an artifact response subtraction process 

9 on the received signals in a computer by estimating the artifact reflection 

10 component of the signal at each antenna location as a combination of the 

11 received signals at the other antenna locations passed through filters to 

12 provide corrected signal data, the filters having weights chosen to 

13 minimize the received signal over that portion of the received signal 

14 dominated by the artifact response. 

1 28. The method of Claim 27 wherein the step of 

2 transmitting microwave signals comprises transmitting microwave pulses 

3 having pulse widths on the order of 100 picoseconds or less in duration. 

1 29, The method of Claim 27 wherein the step of 

2 transmitting microwave signals comprises transmitting microwave pulses 

3 having frequency content at 1 0 GHz or higher. 

1 30. The method of Claim 27 including compensating the 

2 corrected signal data for frequency dependent scattering. 

1 31 . The method of Claim 27 including transmitting the 

2 microwave signals from an array of antennas so as to focus the 

3 microwave power on a candidate location. 

1 32. A method of carrying space-time microwave imaging 

2 of an individual comprising: 

3 (a) transmitting microwave signals from a plurality of 

4 antenna locations into an individual to be examined; 
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(b) 



receiving backscattered microwave signals at a 



6 plurality of antenna locations to provide received signals from the plurality 

7 of antenna locations; 



9 on the received signals in the computer by estimating the artifact 

10 reflection component of the signal at each antenna location as a 

1 1 combination of the received signals at the other antenna locations passed 

12 through filters to provide, corrected signal data, the filter weights chosen 

13 to minimize the received signal over that portion of the received signal 

14 dominated by the artifact response; 

15 id) processing the corrected signal data in a computer in a 

16 beamformer process by time shifting the received signal data to 

17 approximately align the returns from a scatterer at a candidate location, 

18 passing the time aligned signals through a bank of filters with one filter in 

19 each channel^ summing the outputs of the filters to provide a summed 

20 signal, and time-gating the summed signal to provide a time-gated signal, 

21 and calculating the power in the time-gated signal to produce a 

22 beamformer output signal, wherein the filters are designed to present 

23 maximum gain to scattered signals originating from the candidate 

24 location; and 

25 (e) then scanning the beamformer process to a plurality of 

26 different candidate locations in the individual and repeating steps (a), (b) 

27 (c) and (d) at each candidate location by changing the time shifts, filter 

28 weights and time gates for each candidate location to generate multi- 

29 dimensional output power data/ 

1 33 • The method of Claim 32 wherein the step of 

2 transmitting microwave signals comprises transmitting microwave pulses 

3 having pulse widths on the order of TOO picoseconds or less in duration. 



8 



{C> 



carrying out an artifact response subtraction process 
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1 34. The method of Claim 32 wherein the step of 

2 transmitting microwave signals comprises transmitting microwave pulses 

3 having frequency content at 10 GHz or higher. 

1 35. The method of Ciaim 32 wherein the step of filtering is 

2 carried out in the time domain. 

1 36- The method of Claim 32 wherein the step of filtering is 

2 carried out in the frequency domain. 

1 37- The method of Claim 36 wherein the step of filtering 

2 in the frequency domain is carried out for each signal from an antenna 

3 location in a channel by performing a fast Fourier transform on the signal 

4 data to provide frequency domain data, applying weights to the frequency 

5 domain data in each channel and adding the weighted data; and then 

6 performing an inverse fast Fourier transform on the weighted data to 

7 provide a fiitered time domain data signal. 

1 38. The method of Claim 32 including compensating the 

2 received signal data for frequency dependent scattering before carrying 

3 out beamformer processing. 

1 39. The method of Claim 32 including transmitting the 

2 ■ microwave signals from an array of antennas so as to focus the 

3 microwave power on a candidate location. 

1 40. The method of Claim 32 wherein microwave signals 

2 are provided to one antenna at a time and backscattered microwave 

3 signals are received from one antenna at a time for each of the antenna 

4 locations. 
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1 . 41 , The method of Claim 32 wherein microwave signals 

2 are transmitted from ail of the antennas simuJtaneousiy and backscattered 

3 microwave signals are received from all of the antennas simultaneously. 

1 42. The method of Claim 32 including in the beamformer 

2 process the steps of applying a selected window to the time aligned 

3 signals before filtering the signals and further applying a selected window 

4 to the summed signal before calculating the power in the summed signal. 

1 43. The method of Claim 32 wherein the power is 

2 calculated at each candidate location as the sum of the squares of each 

3 sample of the summed signal. 

1 44. The method of Claim 32 wherein filtering is carried out 

2 utilizing FIR filters designed to satisfy a penalized least squares condition 

3 to present unit gain to scattered signals originating from a candidate 

4 location. 

1 45 p The method of Claim 32 wherein the step of 

2 transmitting the microwave signals is carried out simultaneously from all 

3 of the antenna locations by passing microwave pulses for each antenna at 

4 an antenna location through a delay and a filter for each antenna, the 

5 delays and filters selected to focus the radiated microwave energy from 

6 the antennas at a selected candidate location in the object. 

1 46. A method of carrying out hyperthermia treatment of 

2 ah individual with microwave radiation comprising: 

3 (a) providing microwave pulse signals of a selected width 

4 and repetition rate; 

5 (b) providing an array of antennas for radiating 

6 microwave; and 
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7 (c) passing the pulse signals through a separate channel 

8 for each antenna in the array, in each channel selectively delaying the 

9 pulses and filtering the pulses before applying the delayed and filtered 

10 pulses to the antenna for that channel, the delays and the filters selected 

1 1 to focus the radiated microwave energy from the array of antennas at a 

12 selected location m the individual. 

1 47. The method of Claim 46 wherein the pulses have 

2 pulse widths on the order of 1 00 picoseconds or iess in duration. 

1 48- The method of Claim 46 wherein the microwave 

2 pulses have frequency content at 1 0 GHz or higher, 

1 49. An artifact reflection compensation system for use in 

2 space-time microwave imaging comprising: 

3 (a) an array of antennas for radiating and receiving 

4 microwaves; 

5 (b) a microwave pulse source connected to the 

6 array of antennas to provide pulse signals of a selected width and 
. 7 repetition rate to the antennas; 

8 {c) a receiver connected to the antennas to detect 

9 the microwave signals received by the antennas and provide signal data 

10 corresponding thereto; and 

11 (d) a computer connected to receive the signal 

12 data, the computer programmed to estimate an artifact reflection 

13 component of the signal at each antenna as a filtered combination of the 

14 signals at all other antennas and to subtract the estimated artifact 

15 reflection component from the signal data to provide corrected signal 

16 data. 



-61- 



03003907A2 I > 



wo 03/003907 PCT/US02/21397 

1 50. The system of Claim 49 wherein the microwave pulse 

2 source provides pulses having pulse widths on the order of 100 

3 picosecorids or less. 

1 51 - A microwave system for carrying out hyperthermia 

2 treatment of an individual comprising: 

3 (a) an array of antennas for radiating and receiving 

4 microwaves; 

s (b) a microwave pulse source connected to the 

6 array of antennas to provide pulse signals of a selected width and 

7 repetition rate to the antennas; and 

8 (c) signal processing circuitry having a separate 

9 channel for each antenna in the array, the circuitry in each channel 

10 selectively delaying the pulses in a delay and filtering the pulses in a filter 

11 before applying the delayed and filtered pulses to the antenna for that 

12 channel, the delays and the filters selected to focus the radiated 

13 microwave energy from the array of antennas at a selected location in the 

14 individuaL 

1 52- The system of Claim 51 wherein the pulses have pulse 

2 widths on the order of 100 picoseconds or less in duration. 

1 53. The system of Ciaim 52 wherein the microwave pulses 

2 have frequency content at 10 GHz or higher. 
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pulses from a clock (23). The output of the signal generator (21), which may be a short broadband pulses or a signal synthesized 
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each of ihe antennas (35) which are arranged in an array (36) antennas. 
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